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Abstract 

In this report we address the question whether aging in the non equilibrium glassy state is 
controlled by the equilibrium a-relaxation process which occur at temperatures above T g . Re- 
cently Lunkenheimer et. al. [Phys. Rev. Lett. 95, 055702 (2005)] proposed a model for the 
glassy aging data of dielectric relaxation using a modified Kohlrausch-Williams- Watts (KWW) form 
exp — (^age/^~age)^ agc • The aging time tage dependence of the relaxation time T a g e is defined by these 
authors through a functional relation involving the corresponding frequency f(£ age ) = l/(2vrr ag e), 
but the stretching exponent /3 age is same as the /3 a , the a-relaxation stretching exponent. We 
present here an alternative functional form for T ag0 (i age ) directly involving the relaxation time it- 
self. The proposed model fits the data of Lunkenheimer el. al. perfectly with a stretching exponent 
Page different from (3 a . 







Understanding the dynamics of the supercooled liquid in the non equilibrium state has 
been one of the most challenging problems in condensed matter physics. When cooled fast 
enough, the supercooled liquid remains trapped in a specific part of the phase space of the 
constituent particles and cannot equilibrate. For such a system the relaxation time to equi- 
librium increases far beyond the time scale of the experiment. Analysis of the dynamics 
of the liquid in the non equilibrium state reveals a variety of phenomena like aging and 
memory effects 

HQ. 

Theoretical approaches for understanding the complex relaxation 
behavior in the non equilibrium state include computer simulations [3( and study of simple 
dynamical models 0]. An important question in this regard is whether the dynamics in the 
non equilibrium state can be understood as an extrapolation of the alpha relaxation process 
characteristics of the equilibrium states at temperatures above the calorimetric glass tran- 
sition temperature T g I n a recent letter Lunkenheimer et. aZ.[5J have studied the time 



dependent dielectric loss data |6| for various glass formers below the glass transition tem- 
perature. By modeling the experimental data with a modified Kohlrausch- Williams- Watts 
(KWW) form, exp — (t a ge/T a ge(£age)) these authors demonstrate that the non equilib- 
rium dynamics is fully determined by the relaxation times and the stretching parameters 
of the equilibrium a-relaxation process at T > T g . In the present report we demonstrate 
that the data of Ref. [5] can also be fitted in an alternative scheme with a modified KWW 
form same as that proposed in Ref. |5j, but with a relaxation time r age whose time evolution 
is entirely different. Furthermore the proposed model which fits the data of Ref. j|3] per- 
fectly obtains a stretching exponent different from that of the corresponding a-relaxation. 
It therefore conforms to a scenario in which the nature of the relaxation and heterogeneity 
controlling the aging process is different from that of the equilibrium structural relaxation. 
First we consider here the data of dielectric loss e" 5j for glycerol at 179K ( having 

on 

calorimetric glass transition temperature T g = 185K and fragility m = 53 [3, |8j) over the 
frequency range 1Hz - 10 5 Hz. During aging, e" for each frequency decreases continuously with 
t age approaching the equilibrium value over the longest time scale of 10 6 sec. If the aging data 
is fitted with the simple stretched exponential form, we obtain stretching exponent (3 and 
the relaxation time r which are very different from those corresponding to the equilibrium 
a-relaxation at T > T g [9]. The data for the dielectric relaxation function is fitted in this 
case to the form 
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e" Oage) = (4 - 4,) eX P [-(WMge)^ 6 ] + 4j (!) 

where the subscripts "st" and "eq" respectively refer to the initial (t age — > 0) and the long 
time (t age — * oo) limiting values of e" . In such a scheme generally both e'^. and e" eq are 
obtained as the best value for the corresponding fitting parameters. The relaxation time 
T agc and stretching exponent /3 age of the KWW form are also used as free fit parameters. 
The quality of this fitting with the dielectric data from different samples are available in 



Ref. [lCj ( see fig. 1 and 2 there ). Both r age and /3 age as obtained from the fitting of the 
relaxation data in the aging regime show strong frequency dependence. If on the other hand 
we adopt a fitting scheme in which /3 age is kept fixed at its corresponding a-relaxation value 
(3 a = .55, then the data do not fit the form (TjQ) with only r age as a free fit parameter. The 
poorness of such a fitting procedure is discussed further below (see fig. 4). 

An interesting interpretation of this data came from Lunkenheimer et. al. by using 
a modified KWW form with the relaxation time r age being dependent on aging time t age 
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121 ] . However in Ref. 5j the time dependence of r age is prescribed not in terms of the 



relaxation time but the corresponding frequency z/ age which is defined as, 
Lunkenheimer et. al. choose the aging time dependence of v in (T5]) in the following form 



^age (*age) = Ost ~ ^eq) eXp - (27TZ/t age ) /3a ] + U eq . (3) 

According to ([3]), the relaxation time r agc — > l/(27rz/ st ) and l/(27rz/ cq ) as t age — > and oo 
respectively. Lunkenheimer et. al. obtain almost a perfect fit for the dielectric data over 
the whole frequency range with the above choice of the time dependence for the relaxation 
time. They observe that an important feature of this aging process is that the stretching 
exponent /? age is same as that of the a- relaxation f3 a . It therefore implies that the stretching 
of the relaxation remains unaffected by aging below T g and hence conform to the validity 
of the time temperature superposition during the aging process. Furthermore the best fit 
value obtained for u eq in corresponds to a time r eq = l/(27rz/ e q) which agrees with the 
a-relaxation time r a for glycerol extrapolated from higher temperatures ( T > T g ) to sub-T 9 
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regions. This matching of r eq with the extrapolated equilibrium ct-relaxation time values in 
case of glycerol will be further discussed below with fig. 5. 

While this is a very instructive way of interpreting the evolution of the non equilibrium 
state, the anasatz Q for determining the time evolution of r age in terms of the corresponding 
frequency is somewhat unnatural. Instead of the relaxation time itself being time dependent, 
the corresponding frequency v becomes time dependent in ([3]). The latter is then transmitted 
to r ag0 through the standard relation ([2]). As an alternative to this, in the present paper 
we adopt a more natural scheme for the time evolution of the relaxation time r age in the 
following manner. 

Tage^age) = (r st - T fn ) /(t age ) + T fn . (4) 

where the function f(t) reduces to the value 1 or 0, for t — > and oo respectively so that 
the relaxation time r agc attains the asymptotic values r st and Tf n respectively in the above 
two limits. A perfect fit is obtained for all the e (t age ) data at different frequencies using Eq. 
(Cp) with T agc (t age ) is being determined with the ansatz (Tj0). In doing this the choice of the 
function f(t) is not unique. However there are some characteristics that can be associated 
with this function. As indicated above it must change from 1 to as the time changes 
from t = to oo limit. It is plausible to expect the time dependence of f(t) is changing 
monotonically from 1 to 0. We have tested in terms of a reduced time x = t/r, where r 
is the relaxation time the following three options: (a) exponential ( simple or stretched) 
relaxation like exp(— ar) (b) power law type decays x~ a , or (c) a function which decays with 
the tanh(x). Among these the last form gives the best fit to the aging data. We obtain for 
the normalized function f(t), 

f(t) = ^ -„ (5) 

where a Q = 2 1/l3 is a normalization constant to ensure that f(t) reduces to the above stated 
values in the two limits. Note that for times large compared to the relaxation times this 
function also behaves like a stretched exponential function with exponent (3. Therefore the 
function defined in the self-consistent eqn. (jSJ) also represents stretching of the relaxation 
process. The exponent j3 for the stretching process is assumed to be same as the stretching 
exponent /3 age used in the modified KWW form to keep the number of fitting parameters to 
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a minimum. Thus in spite of the differences in the defining relations for the respective r's, 
the /3 age represents stretching of the relaxation process in both the cases (Ref. ^] and ours). 
The fittings (corresponding to different uj's) over the whole frequency range of the dielectric 
data are shown in fig. f . In the present fitting scheme unlike that of Ref. [5| we obtain a 
stretching exponent /3 age = .29 which is different from the corresponding a-relaxation value 
P a = .55. The best fit values for the parameters e" st and ( at a given frequency ) obtained 
respectively in the present work and Ref. [5J agree closely. This comparison is displayed 
in fig 2, We display next in fig. 3 the aging time dependence of r age in our calculation and 
that of Ref. j^j. In the present case the the relaxation time is strongly time dependent 
at the initial stage and saturates at relatively earlier as t age reach the time Tf n (say). In 
our fitting scheme the best fit result value for Tf n is very different from the corresponding 
T cq of Lunkenheimer et. al. and cannot be simply related to the a relaxation process. 
All the dielectric data for e scale into a single master curve as shown in the fig. 4. In 
the same figure we display the curve corresponding to the fitting of Lunkenheimer et. al. 
being practically indistinguishable from ours. We also display here for comparison the best 
fittings curves which are obtained with simple KWW function having constant relaxation 
times, respectively equal to r st and r eq . The stretching parameter in both cases is (3 a . The 
inadequacy of a simple KWW form in fitting the observed data is clearly observable here. 

In addition to the dielectric data for glycerol, following Lunkenheimer et. al. we also fit 
using our scheme the dielectric data for Propylene carbonate ( PC, T g fh I59K, m = 104 
8]), Propylene glycol( PG, T g « 168K, m = 52 3), xylitol ( T g « 248K, m = 94 Q,Q) 
and structural relaxation for [Ca(NOs)2]o.4[KN03]o.6 ( CKN, T g ps 333K, m = 93 ) and the 
results are respectively shown in fig 5a-5d. The temperature ( below the corresponding T g 
) in each case is also displayed in the figure. The variation of the relaxation time T age (t age ) 
in the modified KWW form for fitting the data in each of the above materials is shown in 
fig 6. For all these materials, the general trend in aging time dependence of the relaxation 
time is the same as that in case of glycerol. r ag0 is strongly time dependent initially ( up 
to time Tf n ) and finally saturate to an almost constant value as is clearly visible in fig. 6. 
The TABLE I represents the stretching exponent /3 age as a best fit parameter in the aging 
data for different materials obtained from present work (column 1) and also from Ref. 5| 
(column 2). 

In case of the fitting scheme of Lunkenheimer et. al. the time scale r eq agree with the 
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extrapolated value of the a-relaxation to the sub T g region. This is shown in fig. 7 for all 
the five materials whose relaxation data have been considered above. First r^'s for these 
materials with inverse temperature 1/T are shown. On the same figure we display for each 
curve ( with a filled circle ) the corresponding value of r eq = l/(2m/ eq ) obtained from the best 
fit value of ([3]) obtained by Lunkenheimer et. al. These points all lie on the corresponding 
extrapolated curve of the a-relaxation times in the sub-T 9 regime. These authors conclude 
that the equilibrium a-relaxation is determining the aging process . On the other hand the 
corresponding values of r eq = Tf n obtained in our fitting scheme are shown for each curve 
with a filled triangle conforming to a relaxation mechanism of aging different from that of 
the equilibrium a-relaxation. 

We have presented here an alternative scenario for explaining the dielectric relaxation 
data of Lunkenheimer et. al. [s] in the aging regime. The present model proposes that the 
aging process involve two basic steps. In the first stage, the aging data fits to a modified 
stretched exponential form with a time dependent relaxation time r age . This is similar to 
the scheme of Lunkenheimer et. al. but the time dependence of T age (t age ) is more natural 
here. As t age reach a characteristic time scale Tf n the data can now be described in this 
second stage with a simple KWW form having a constant relaxation time comparable to 
Tf n . This holds simultaneously for relaxation data at all frequencies as shown in fig 6. The 
limiting value of Tf n obtained here from fitting the data is however not same as that obtained 
from extrapolation of a- relaxation times at higher temperatures to the sub-T 5 region. The 
stretching exponent /3 age is also very different from that of the a-relaxation process. The 
time dependence of the aging process and the corresponding relaxation time is possibly 
controlled by mechanisms which are different from that of equilibrium a-relaxation. 

The present analysis do not to claim in any way that the fitting scheme proposed by 
Lunkenheimer et. al. is invalid. The main difference of the model proposed here from that 

n 

in ref. [5|] lies in the manner in which the relaxation time T age depends on t age . In our scheme 
the relaxation time r ag0 is directly dependent on £ age . This is distinct from adopting a some 
what unusual scheme of ([3]) in which the time dependence is imposed in the frequency u age . 
The justification of this scheme in which aging time dependence of r age is defined in terms 
of the frequency is not obvious. The motivation for adopting such a scheme is that the 
asymptotic value for r age (for long times) agrees with value of a-relaxation times at higher 
temperatures extrapolated to the sub-T 9 region. There is no compelling reason to assume 
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such a link between the non-equilibrium aging process and the equilibrium a-relaxation. 

To explain the last statement let us consider the correlation function C(t + t age ,£ ag0 ) of 
fluctuations at two different times t age and t + t age . In the non-equilibrium state C depends 
on two times, i.e., both t and t age . When time translational invariance holds, the dependence 
of C on t age disappears. This happens over some time scale which is r eq . The correlation 
function now depends only on t and the system has reached equilibrium. At this stage 
the relaxation is controlled by the a- relaxation time r a . But the time scale r eq which sets 
the dependence of the correlation function on t age has no reason to be same as r a . In other 
words we insist that there is no aprirori reason to assume that the time evolution of e during 
the aging process in the non equilibrium state is controlled by the equilibrium a-relaxation 
process and the corresponding heterogeneity. Our work clearly shows here that the dielectric 
data can also be fitted with the modified KWW form with a different choice for the aging 
time dependence of the corresponding relaxation time. 

An important aspect of the proposed relaxation behavior lies in the different time evolu- 
tion of the relaxation time r in comparison to that of Ref. 5J . Here r decreases with waiting 
time t age , implying that the aging process accelerates with time. It should however be noted 
that in the present case the relaxation time only decrease initially and then becomes almost 
constant. The stretching data at different frequencies correspond to a single relaxation time 
which is ahnost constant for aging time ^ beyond this initia. sca.e r,„. The stretching 
exponent /3 age for different frequency data is also the same. In Ref. [5|] the relaxation time r 
actually grows with waiting time. Though somewhat speculative at this stage, such a differ- 
ence in the basic nature of the aging process might imply a different microscopic mechanism 
for aging altogether. The more appropriate choice between the two schemes ( for fitting the 
aging data) discussed here can only be ascertained with a proper theoretical model for the 
evolution of the non equilibrium state. 
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Materials 


Stretching exponent /3 age 


Present Work 


Ref. [5] 


Glycerol 




55 f5l 


PC 


0.24 


0.6 [6,15] 


PG 


0.23 


0.58 [16,rT] 


xylitol 


0.45 


0.43 [16] 


CKN 


0.34 


0.4 [18] 



TABLE I: Tne Aching exponent a, . best * palter in tne a gmg data f„ r different m a- 
terials obtained (a) from present work (column 1) (b) from Ref. [5j (column 2). The corresponding 
reference from which the data for a-relaxation is taken are also provided in col. 2 
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FIG. 1: e (t a ge) data for glycerol at 179°K vs. i age (in sec) at the frequencies displayed in the inset. 
Solid lines correspond to the best fits using the modified KWW form with T age (£age) determined 
from eqn. (U]). 
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FIG. 2: The best fit values of e st for the curves in fig. 1 vs. frequency zv(Hz) as obtained with the 
present work (solid line) and the corresponding results of Ref. 5] (dashed). The inset shows the 
best fit e" with frequency v respectively from the present work and Ref. 5(. 
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FIG. 3: The best fit r(t age ) vs. i age (in sec) for glycerol using ansatz Q of the present work. 
Result of Ref. [5[ for r(t aS e) as inset. 
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FIG. 4: Scaling of the dielectric data for glycerol at 179°K corresponding to different frequencies. 
Fits obtained with T age (i age ) in the modified KWW ansatz determined from (a) Eqn. (|4]) of the 
present work (solid line) and (b) from Ref. (dashed line). Also shown are best fits obtained 
taking r age = r s t (short dashed) and r age = r eq (long dashed) respectively with the stretching 
exponent being fixed at f3 a . 
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FIG. 5a: e" (t age ) data for PC 
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at 152° K vs. t age (in sec) at the frequencies displayed in the 



inset. The solid lines are obtained using our scheme. 
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FIG. 5b: e"(t age ) data for PG 8 
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16l | at 157° K vs. t age (in sec) at the frequencies displayed in the 



inset. The solid lines are obtained using our scheme. 
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FIG. 5c: e (t age ) data for xylitol [iff] at 243° K vs. i age (in sec) at the frequencies displayed in the 
inset. The solid lines are obtained using our scheme. 
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FIG. 5d: M" (i age ) data for CKN 



181 ] at 329°K vs. t agc (in sec) at the frequencies displayed in 



the inset. The solid lines are obtained using our scheme. 
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FIG. 6: r(f age ) vs. i age (in sec) for different materials using ansatz in the present work. 
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FKj jj: , m e r . v, inverse of „ Ure for CKN H , xvlitolQ, glvceroQ, 

PG[17], PC 18]. The T e q obtained from the fitting of the corresponding data is shown by filled 
circles for the fitting of Ref. [5j and by filled triangles the present fitting scheme. 
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